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We propose an entanglement mechanism of nuclear spins in quantum dots driven by the 
electric current accompanied by the spin flip. This situation is relevant to a leakage current 
in spin-blocked regions where electrons cannot be transported unless their spins are flipped. 
The current gradually increases the components of larger total spin of nuclei. This correlation 
among the nuclear spins markedly enhances the spin-flip rate of electrons and hence the leakage 
current. The enhancement of the current is observable when the residence time of electrons in 
the quantum dots is shorter than the dephasing time T2 of nuclear spins. 

KEYWORDS: quantum dot, entanglement, nuclear spins, hyperfine interaction, spin blocl<ade, dephasing, super- 
radiance, Dicke effect 
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The spin relaxation of electrons in quantum dots is 
an important issue. The relaxation time must be sufR- 
ciently long for the implementation of quantum comput- 
ing devices utilizing electron spins in quantum dots.^^ 
The time has been examined experimentally by opti- 
cal^' and transport measurements.'*'^-' Mechanisms of 
the spin relaxation have been investigated theoretically 
by, 6.17., hyperfine interaction,^"*^) spin-orbit interac- 
tion*^' *'^) and higher-order tunneling processes (cotun- 
neling). 14.15) 

In this letter, we study the unique properties of nu- 
clear spin states in quantum dots, which stem from the 
spin relaxation of electrons through hyperfine interac- 
tion. In quantum dots, an electron occupying a partic- 
ular level couples with numerous nuclear spins simulta- 
neously through the interaction. Hence a spin flip of the 
electron results in a correlation among the nuclear spins. 
The correlation is an entanglement: we do not know 
which of the nuclear spins is flipped. We show that (i) the 
electric current accompanied by the spin flip significantly 
increases the correlation among nuclear spins, and (ii) 
the correlated state of nuclear spins markedly enhances 
the spin-flip rate of electrons. The direct dipole-dipole 
interaction between the nuclear spins is much smaller 
than the hyperfine interaction^-* and is thus disregarded. 
Its influence is discussed later. 

This situation is relevant to a leakage current in the 
spin-blocked region where electrons cannot be trans- 
ported unless the spins are flipped. In weakly coupled 
double quantum dots, Ono et al. have observed a cur- 
rent suppression when two electrons occupy the lowest 
energy level in each dot with parallel spins an electron 
tunneling from one dot to the other is forbidden by the 
Pauli exclusion principle (see Fig. 4). In this region, we 
predict that a leakage current accompanied by the spin 
flip in one of the dots is signiflcantly enhanced by the 
nuclear spin entanglement. If our theory is verifled, it 
would be the flrst observation of entanglement in solid- 
state devices. 

Our mechanism is analogous to the current-induced 
dynamic nuclear polarization (DNP) in quantum Hall 

systems. *^"*^) The DNP is created by electron scattering 
between spin-polarized edge states, accompanied by the 
spin flip. The DNP influences the electronic state and, 
in consequence, leads to a hysteresis in the longitudinal 
resistance.^"' Note that nuclear spins are entangled 
but not polarized in a specific direction in our case. 

We examine a simple model for the spin-blocked quan- 
tum dots and discuss its relevance to a realistic situation 

* E-mail address: eto@rk.phys.keio.ac.jp 



later. Our model is as follows. A quantum dot is con- 
nected to external leads, L, R, through tunnel barriers. 
The Coulomb blockade restricts the number of electrons 
in the dot to be A'ci or iVci -f- 1. N^i electrons form a 
background of a spin singlet. From lead L on the source 
side, an extra electron tunnels into the dot and occu- 
pies a single level with an envelope wavefunction ili{r). 
The spin of the electron is either | t) {Sz = 1/2) or 
I i) {^z = ~l/2) with equal probability (we assume an 
easy-axis of electron spins in the z direction, e.g., in a 
magnetic field). The electron stays for a long time by 
the spin blockade, interacting with N nuclear spins, Ik 
(located at r^; k = I,-- - ,iV), through the hyperfine 
contact interaction. After the spin is flipped, the elec- 
tron tunnels out to lead R on the drain side, and the 
next electron is immediately injected from lead L. 
The hyperfine interaction is expressed as 



Hm — 2 



N 

E 

fc=i 



akS ■ Ik, 



(1) 



where ak oc |V'('''fc)P- ^ ™ GaAs quantum dots. 
We consider nuclear spins of 1/2 for simplicity. We as- 
sume that the spin-fiip rate from | J,) to | f) (| f) to | J,)) 
is given by the perturbation of i/hf , 



r-A|(vi/fi„;ta)|i/M|vi/init;i (?)>!' 



(2) 



where 5'init and are the initial and final states of 
nuclear spins, respectively.^^-' We first study the case of 
homogeneous coupling constants, ak = ex.. This is a good 
approximation for the majority of nuclear spins since the 
distance between neighboring nuclei is much smaller than 
the size of the quantum dot. Then the Hamiltonian (1) 
is reduced to 



(3) 



where I = X^^i -^fc ^he total spin of nuclei. A generic 
case with inhomogeneous hyperfine coupling is examined 
later. 

The Hamiltonian (3) [and also (1)] indicates that N 
nuclear spins interact with a common "field" of an elec- 
tron spin although there is no direct interaction between 
them (dipole-dipole interaction between nuclear spins is 
weak and neglected). This field results in an entangle- 
ment among nuclear spins. To illustrate this, let us con- 
sider the simplest case of = 2 and begin with a po- 
larized state of nuclear spins \Iz = l/2)i\Iz — 1/2)2- 
An electron with spin | J,) tunnels into the dot and is 
spin-fiipped. Then the state of nuclear spins becomes 
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(I - l/2)i|l/2)2 + |l/2)i| - 1/2)2)/a/2 ^ \J=l,M^ 0), 
where J and M are the total spin and its z component, 
respectively. This is an entangled state; we do not know 
which of the nuclear spins is flipped. After the electron 
tunnels out of the dot, the next electron is injected with 
I J.) (or I t)) and interacts with this state of nuclear 
spins. The spin-flip rate for the second electron is given 
by 2Aa^ . This value is twice the rate in the case of 
non-entangled states of nuclear spins, |l/2)i| — 1/2)2 or 
I — 1/2)1 1 1/2)2. In general, the capability of state \ J,M) 
to flip an electron spin is Aa'^{J ± M)(J =f M + 1). 

Now we formulate the case of N nuclear spins. We 
assume that the initial state of the nuclear spins is 



= 



i=±l/2 /;j,jv=±l/2 



x\Iz,i}i\Iz,2)2---\Iz,n} 



N, 



(4) 



where {C{Iz.i--'' - Izm)} f^rc randomly distributed.^'^) 
By the transformation of the basis set into the cigcn- 
states of the total spin J and its z component M of all 
the nuclei (J < N/2, - J < M < J), eq. (4) is rewritten 
as 

^(0)^ E C^°ljJ^M,X), (5) 

J,M,A 

with random coeflicients {c';"ljj^} (E,/,m,a I'^'j.mJ^ = 
1). Index A distinguishes states with the same J and M. 
The number of such states is 

K{N, J) = (2J + l) ^ ^+2^+2^ ;(iV-2j y- (6) 

The first electron with | J,) (or | f)) enters the dot and 
interacts with ^f'"' of nuclear spins. The spin-flip rate is 
= Aa^i?^^, where 

^(0)^ E \C^]!lj\J±M){JTM+l). (7) 

J,M,X 

When > 1, \cf^it^^\^ can be replaced by 1/2^ (law of 

large number; fluctuation of 1/V2^ is neglected). Then 
= (1/2^)E,/,M W^)(^±M)(JtM + 1) « 

N/2. Note that T^") = {N/2)Aa'^ is i dentical to the spin- 
flip rate which would be evaluated under the assumption 
that one of the nuclear spins is flipped with the electron 
spin.24) 

After the spin flip, the electron tunnels ofi^ the dot. 

The state of nuclear spins becomes 



1 



x|J,Mt1,A). (8) 



^l'^^' includes more components of larger J and smaller 
\M\. The correlation among the nuclear spins increases 
each time an electron is injected, spin-flipped, and 
ejected out of the dot. After n such events, the state 
of nuclear spins is 



.J,M.X 



cfi,|J,M,A), 



(9) 



{J±M){JtM + 1) 



■C^Hm'I (10) 
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Fig. 1 . Distribution of (a) total spin J and (b) its z component in 
the state of nuclear spins, (a) p{J) = A I*-^j"m aI^' ^^'^ (^) 
p{M) = X)j A l^j^M aI^' where n is the number of transported 
electrons accompanied by the spin flip. The number of nuclear 
spins is N = 100. In eq. (10), we taJie the geometrical mean 
between upper and lower signs. 
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in) 



^ l4i,|^(J±M)(JTM + ^).l) 



J,M,A 



is expressed as = /(") //(""^^ with 
/^"^ = 4 E ^(^' [(-^ ±M){JtM + ir . (12) 

J,M 

Figure 1 shows the distribution of (a) the total spin J 
and (b) its z component M. With increasing n, the 
weight of larger J and smaller \M\ increases. This means 
that the total nuclear spins are developed in the plane 
perpendicular to the easy-ax;is of electron spins. 

The correlation among the nuclear spins enhances the 
spin-flip rate. The rate for the (n+ l)th electron is given 
by r(") = Using eq. (12), we find that « 

(A/2)n for K n < A/2, and F(") « {N/2)^ for A/2 < 
n. Therefore the spin-flip rate increases with n linearly 
(r(") « nr(0)) and finally saturates (F^") w (A/2)r(")j, 
as shown in Fig. 2(a). 

We evaluate the electric current I{t). The time interval 
between nth and (n-l- l)th electron transports is At^") 
l/r("). The current as a function of t, I{t) = eF^") at 
t = X]"=o is shown in Fig. 2(b). The approximate 

form of F*^") yields an asymptotic form of the current, 



I(t) 

m 



er(o)er""* 
e(A/2)F(o) 



{t < tsat) 

(isat <i)- 



(13) 



The current grows significantly with time and finally 
saturates. The saturation time is given by fgat = 
ln(A/2)/F(°), where F(°) is the spin-flip rate of an elec- 
tron with the uncorrected state of nuclear spins. 

Up to now, we have examined the case of homoge- 
neous hyperfine coupling in the quantum dot, ctfc = a. 
Generally, the coupling constant ak depends on the po- 
sition of the nucleus due to the spatial variation of i^{r). 
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Fig. 2. (a) Spin-flip rate F'^"' as a function of n (number of trans- 
ported electrons accompanied by spin flip) and (b) electric cur- 
rent I{t) as a function of time t, with homogeneous hyperfinc 
coupling. The number of nuclear spins is Af = 100. Wc take the 
geometrical mean between upper and lower signs in eq. (12). 
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Fig. 3. Spin-flip rate F^"' as a function of n (number of trans- 
ported electrons accompanied by spin flip), when the hyperfine 
coupling constants are aj. oc cxp [— (fe — X)lxo\ (fc = 1, 2, ■ ■ ■ , N\ 
i J] J, = q). A'' = 14 with 2:0 = 10 (solid lines) or 00 (broken 
line). The residence time of the first electron in the dot, t, is 

(a) oltIK = 0, (b) 10, (c) 20, and (d) 40. F(") is averaged over 
the random distribution of initial states. The error bars denote 
the variation of F'"^. 



energy conservation in the spin flip.^' The Hamiltonian 
in the dots is H = + Hp^ + Ht + Hm + ifei-ph, where 



With this inhomogeneous couphng, (i) the total spin J 
is no longer a good quantum number, and (ii) the cor- 
related state of nuclear spins is dephased. During the 
stay of an electron in the dot [r ~ l/r(°) = 2/{Aa^N)\, 
each term in eq. (4) gains a phase factor e"^"^ with 
w = ±y2k^klz,k/f>'- This factor diminishes the corre- 
lation ifthe residence time t is sufficiently long. 

To examine those cifFects. wc perform numerical cal- 
culations in the case of at oc exp— [{k — l)/xo]^ {k = 
1, 2, • • • ,N; jf^f.cxk = a). The number of nuclear spins 
is W = 14. Electrons with spin | f) and | |) are injected 
into the dot alternately. For the initial state of nuclear 
spins, the random average is taken over the coeffi- 

cients {C{Iz,i, ■ ■ ■ ,Iz,n)} in eq. (4). 

Figure 3 presents the spin-flip rate F*^"^ with a;o = 00 
(homogeneous coupling; broken line) and with xq = 10 
and ar/h = 0, 10, 20, and 40 (solid lines). When the 
dephasing effect is negligible (r = 0), the spin- flip rate 
is significantly enhanced by the entanglement of nuclear 
spins, even with finite Xq- The behavior of F^") is char- 
acterized by the effective number of nuclear spins which 
couple to the electron spin (iV N^g xq in the previ- 
ous discussion). With an increase in the residence time 
of electrons in the dot, r, the enhancement of F^") be- 
comes less prominent, and is negligible when r is larger 
than the dephasing time of nuclear spins by the inhomo- 
geneous hyperfinc field. 

In actual quantum dots, there arc other origins of the 
dephasing of nuclear spins, e.g., dipolc-dipole interac- 
tion. Generally, as long as the dephasing time T2* exceeds 
T, our mechanism appears robust against the dephasing 
effects because the spin flip of electrons drives the nuclear 
spins to be correlated every time. 

Our model is relevant to the experimental situation of 
double quantum dots shown in Fig. 4. A spin blockade 
takes place when the spin of an incident electron in dot 
L is parallel to that of an electron trapped in dot R.^^ 
We assume a magnetic fleld in which the Zeeman energy 
for electrons, Ez, is much larger than a in iJhf, whereas 
the Zeeman energy for nuclear spins is negligible. The 
electron-phonon interaction is taken into account for the 



E 

j=L,R 



E( 



el— ph 



where aj aj^cr {b\, bq) are creation and annihilation op- 
erators, respectively, for an electron in dot j (a phonon). 



a]^aO-jM and 



In the spin blockade of |L T)l^ T)- an electron tun- 
neling to \R i)\R t) is accompanied by the spin flip in 
one of the dots, e.g., dot L, and emission of a phonon of 
hojq = 5E = sl ~ {er. + U) — Ez. The transition rate is 



F = A|(4t';il^hf|*L1(;T)Pwith 



^ ^ 2^ /y /3L,q - 0R.. 
h V EzSE 



DphiSE) [NphiSE) + 1] 



in the lowest order of Ht + iJhf + i?ei-ph- Here, Dph is 
the density of states for phonons and A^ph is the number 

of phonons. This agrees with eq. (2). 

Besides eq. (2), wc have assumed that the nuclear spin 
state is a pure state after an electron tunnels out from 
dot state \R l)\R t) to lead R and is detected by cur- 
rent measurement. For this assumption to be valid, 
the following higher-order tunneling processes (cotun- 
neling)^' ^'"'^ must be disregarded: (i) third-order pro- 
cesses in which an electron tunnels out from dot R to 
lead R, an electron tunnels from dot L to dot R, and 
finally an electron tunnels into dot L from lead L, or (ii) 
second-order processes in which an electron tunnels out 
from dot L (R) to lead L (R), and another electron tun- 
nels in from lead L (R) to dot L (R). In the presence 
of such processes and hyperfine couplings, the nuclear 
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Fig. 4. An experimental situation of double quantum dots. A spin 
blockade takes place when the spin of an incident electron in 
dot L is parallel to that of the electron trapped in dot R.^^ 
The leakage current in the spin-blocked region is relevant to our 
theory. 



spin state would be entangled with electron states in the 
leads. Then our calculations would no longer be justified. 

We have also neglected the spin-orbit (SO) interaction, 
which may play a role in the spin relaxation of electrons 
in quantum dots.^^' The coexistence of hyperfine and 
SO interactions would complicate the evolution of the 
nuclear spin state although the state remains pure (a lin- 
ear combination of spin-flipped and -unflippcd states) in 
this case. Its analysis is beyond the scope of the present 
study. 

In conclusion, we have proposed a current-induced en- 
tanglement of nuclear spins in a quantum dot. The cur- 
rent accompanied by spin flip in the quantum dot grad- 
ually increases the correlation among nuclear spins when 
the residence time of electrons in the dot does not exceed 
the dephasing time of nuclear spins Tj*. The correlated 
state of nuclear spins significantly enhances the spin-flip 
rate of electrons. A relevant situation to this theory is 
the leakage current in spin-blocked regions. 

The time dependence of the leakage current /(<) shown 
in Fig. 2(b) may be observed in the experiment using a 
sequence of NMR pulses. The nuclear spins are random- 
ized by each pulse. By tuning the interval of the pulses. 
At, /(At) can be measured. Quantitatively, the dephas- 
ing time of nuclear spins is estimated to be Tj* '^-^ 1 ns in 
a case of GaAs quantum dots.*' The saturation time of 
I{t) is of the order of the residence time of electrons in 
quantum dots [tsat — ln(iV/2)T in eq. (13)], which should 
be shorter than T^j*. 

Finally, we comment on an analogy between our mech- 
anism and the Dickc effect of superradiance.^^'^*' In the 
spontaneous emission of photons from N atoms with two 
levels (pseudo-spin Sz = ±1/2), all the atoms could in- 
teract with a common electromagnetic field. The emis- 
sion of photons is significantly enhanced if N atoms are 
excited initially. This is due to the formation of the 
pseudo-spin state |J, M) with J = N/2. Starting from 
I J, J), the state of TV atoms changes like a cascade, | J, J), 
I J, J — 1), • • • , emitting photons rapidly. A similar ef- 
fect has been proposed for the emission of phonons from 
N equivalent quantum dots.^*'^^-* The atoms (quantum 
dots) correspond to the nuclear spins in our model, and 
the emission of photons (phonons) to the spin flip of elec- 
trons. A major difference is the initialization. N excited 
states must be prepared by pumping in the Dicke effect, 
whereas such initialization is not necessary in our mech- 
anism. 

The authors are indebted to L. P. Kouwenhoven for the 
suggestion of the experiment using NMR pulses. They 
gratefully acknowledge discussions with K. Kawamura, 
R. Fukuda, S. Komiyama, T. Inoshita, A. Shimizu, K. 
Ono, G. E. W. Bauer, and S. I. Erlingsson. This work 



was partially supported by a Grant-in- Aid for Scientific 
Research in Priority Areas "Semiconductor Nanospin- 
tronics" (No. 14076216) of the Ministry of Education, 

Culture, Sports, Science and Technology, Japan. 



10 
11 
12 
13 
14 

15 
16 

17 

18 

19 

20 



21 



22 



23) 



24) 

25) 
26) 



27) 
28) 

29) 



D. Loss and D. P. DiVincenzo: Phys. Rev. A 57 (1998) 120. 
M. Paillard, X. Marie, P. Renucci, T. Amand, A. Jbeli and J. 
M. Gerard: Phys. Rev. Lett. 86 (2001) 1634. 
A. S. Lenihan, M. V. Gurudev Dutt and D. G. Steel: Phys. 
Rev. Lett. 88 (2002) 223601. 

T. Fujisawa, D. G. Austing, Y. Tokura, Y. Hirayama and S. 
Tarucha: Nature 419 (2002) 278. 

K. Ono, D. G. Austing, Y. Tokura and S. Tarucha: Science 
297 (2002) 1313. 

S. L Erlingsson, Yu. V. Nazarov and V. I. Fal'ko: Phys. Rev. 
B 64 (2001) 195306. 

S. I. Erlingsson and Yu. V. Nazarov: Phys. Rev. B 66 (2002) 
155327. 

I. A. Merkulov, Al. L. Efros and M. Rosen: Phys. Rev. B 65 
(2002) 205309. 

Y. B. Lyanda-Geller, I. L. Aleiner and B. L. Altshuler: Phys. 
Rev. Lett. 89 (2002) 107602. 

A. V. Khaetskii, D. Loss and L. Glazman: Phys. Rev. Lett 
88 (2002) 186802. 

J. Schliemann, A. V. Khaetskii and D. Loss: Phys. Rev. B 66 
(2002) 245303. 

A. V. Khaetskii and Yu. V. Nazarov: Phys. Rev. B 61 (2000) 
12639. 

A. V. Khaetskii and Yu. V. Nazarov: Phys. Rev. B 64 (2001) 
125316. 

D. V. Averin and Yu. V. Nazarov: Phys. Rev. Lett. 65 (1990) 
2446. 

M. Eto: Jpn. J. Appl. Phys. 40 (2001) 1929. 

B. E. Kane, L. N. Pfoiffer and K. W. West: Phys. Rev. B 46 
(1992) 7264. 

K. R. Wald, L. P. Kouwenhoven, P. L. McEuen, N. C. van 
der Vaart and C. T. Foxon: Phys. Rev. Lett. 73 (1994) 1011. 
D. C. Dixon, K. R. Wald, P. L. McEuen and M. R. Melloch: 
Phys. Rev. B 56 (1997) 4743. 

T. Machida, S. Ishizuka, T. Yamazaki, S. Komiyama, K. Mu- 
raki and Y. Hirayama: Phys. Rev. B 65 (2002) 233304. 
S. Kronmiiller, W. Dietsche, K. von Klitzing, G. Denninger, 
W. Wegscheider and M. Bichler: Phys. Rev. Lett. 82 (1999) 
4070. 

J. H. Smet, R. A. Deutschmann, W. Wegscheider, G. Abstre- 
iter and K. von Klitzing: Phys. Rev. Lett. 86 (2001) 2412. 
We assume that the final state of nuclear spins is a pure state, 
not a mixed state with electron states | T) or | |), after the 
electron tunnels out of the dot. We observe a spin flip of 
electron by the leakage current in a spin-blocked situation, 
whereas we cannot observe which of the nuclear spins was 
flipped. In a spin-flip process, a phonon may be emitted. How- 
ever, if the wavelength of the phonon is much longer than the 
nucleus-nucleus distance, we do not obtain such information 
on the nuclear spins. 

is different from fW' = (ci |l/2)i -|- di | - l/2)i) (g)- ■ ■ ® 
(c]v 1 1/2) jv + df4 1 — 1/2) AT ) , with random coefficients {c^ , } , 
in which the nuclear spins are not entangled at all. Numerical 
calculations using iffO)', however, yield qualitatively the same 
results as those in Fig. 3 obtained using See M. Eto, T. 

Ashiwa and M. Murata: in preparation. 

Hence the entanglement of nuclear spins cannot be observed 
by optical experiments in which the spin relaxation of single 
excitations is measured.^' '^^ 

A similar model has been studied by T. Brandos and B. 
Kramer: Phys. Rev. Lett. 83 (1999) 3021. 
For the transition rate from \L [)\R I) to \R '\)\R J,), &E is 
replaced by SE' = 5E+2Ez - Note that the spin of an incident 
electron in dot L is parallel or anti-parallel to the spin of an 
electron trapped in dot R, with equal probability. In the latter 
case, the spin blockade does not occur; the electron tunnels 
in the second dot and finally out to lead R immediately. The 
current is the sum of the currents with and without spin flip. 
R. H. Dicke: Phys. Rev. 93 (1954) 99. 

T. Brandes: Interacting Electrons in Nanostructures, eds. R. 
Haug and H. Schoeller (Springer, Berlin Heidelberg, 2001). 
T. Brandes, J. Inoue and A. Shimizu: Phys. Rev. Lett. 80 
(1998) 3952. 



